Tuning the exchange interaction by electric field in laterally coupled quantum dots 
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The effect of external electric field on the exchange interaction has been studied by an exact 
diagonalization method for two electrons in laterally coupled quantum dots (QD's). We have per- 
formed a systematic study of several nanodevices that contain two gate-defined QD's with different 
shapes and sizes located between source and drain contacts. The confinement potential is modelled 
by two potential wells with a variable range and softness. In all the considered nanodevices, the 
overall dependence of exchange energy J on electric field F is similar, i.e., for low fields J increases 
with increasing F, for intermediate fields J reaches a maximum, and rapidly falls down to zero if F 
exceeds a certain critical value. However, the J(F) dependence shows characteristic properties that 
depend on the nanodevice geometry. We have found that the low- and intermediate-field behaviour 
can be accurately parametrized by linear function J(F) — aF + f3, where a is independent of the 
nanodevice geometry and softness of the confinement potential. We have shown that the linear 
J(F) dependence appears only if the tunnel coupling between the QD's is weak, i.e., the interdot 
separation is sufficiently large. The J(F) dependence becomes non-linear for the strong interdot 
tunnel coupling. If the QD located near the contact, to which the higher voltage is applied, possesses 
the elliptic shape and is larger than the other QD, the J(F) dependence shows a plateau in a broad 
electric-field regime. The linearity and rapid jumps of J(F) as well as the existence of the plateau 
can be applied to tune the exchange interaction by changing the external electric field. 

PACS numbers: 73.21.La,03.67.Lx 



I. INTRODUCTION 

Exchange interaction is one of the most characteris- 
tic quantum effects in many-electron systems. In natural 
atoms and molecules, it leads to a singlet-triplet split- 
ting of energy levels and is responsible for a binding of 
atoms into a molecule. In solids, it gives rise to a co- 
valent bonding of elemental semiconductors and ferro- 
magnetic properties of metals. In the absence of external 
fields, the exchange interaction in natural atomic systems 
is determined by the charges of nuclei and the number of 
electrons and is essentially fixed. Man-made solid-state 
analogs of atomic systems, namely, quantum dots (QD's), 
called also artificial atoms, and coupled QD's (artificial 
molecules) can be fabricated in various designed shapes 
and sizes. One can also change the depth and range of 
the potential confining the electrons. This gives us a 
unique opportunity of engineering the quantum states of 
electrons confined in the QD's and tuning the exchange 
interaction. The exchange interaction between electrons 
in QD's has been proposed as an effective mechanism for 
changing the electron spin, i.e., for perfo rming the quan- 
tum logic operations on spin qubits This mechanism 
seems to be very promising in quantum information pro- 
cessing with the solid-state nanodevices. 

The investigation of the exchange interaction in the 
QD-based nanodevices allows us to elaborate the meth- 
ods of controlling and tuning this interaction, which in 
turn leads to the controlled manipulation of the elec- 
tron spin qubits f321214| Recently, the exchange-interaction 
induced spin swap operations in coupled QD's have 
been simulated by a direct solution of a time-dependent 
Schrodinger equation^ The coherent manipulation of 



spin qubits in lateral QD's has been studied experimen- 
tally by Petta et al.P Elzerman et al.J^ and Hayashi 
et alP Hatano et alP determined the tunnel and ex- 
change couplings in laterally coupled vertical QD's. The 
quantum logic operations can also be performed in the 
nanowire double QD'sP^ 

The exchange energy is denned as 

J = E T -E Sl (1) 

where Et and Eg are the lowest triplet and singlet energy 
levels, respectively. In physics of solid-state nanodevices, 
exchange energy plays a two-fold role: (i) it deter- 
mines the strength of the exchange interaction (H eisen- 
berg interaction), which can swap the spin qubits EMHI1 
(ii) J can be treated as the exchange splitting (singlet- 
triplet splitting) that allows us to distinguish the differ- 
ent spin states of the electron system. 4 In both the cases, 
we should know how to tune J by applying the external 
fields. 

In the coupled QD system, the ex- 
change energy was calculated by s everal 
groupc . 1 ^^ 1 ' 12 '"' 14 ' 15 ' 16 ' 17 ' 18 ' 19 ' 20 ' 21 ' 22 ' 23 ' 24 ' 2 ^' 2 ' 51 ^ The 

effect of an externa l magnetic field was investigated in 
papers jH 3 l 17 l 18 l 21 l 2 H ^he magnetic field, applied perpen- 
dicular to the plane of the electron movement, lowers 
the energy of the triplet state, which decreases the 
exchange splitting. The asymmetry of the QD's gives 
rise to an enhancement of the exchange interaction! 17 ^ 22 ! 
The si ze effect s in the exchange coupling were studied in 
papers {I 4 I 25 I 27 I The increasing size of the coupled- QDsys- 
tem leads to the decrease of the exchange energy E l 25 l 27 l 
For the two identical elliptic QD's Zhang et al ! 22 1 23 1 
calculated the exchange energy as a function of aspect 
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ratio r — R y /R x , where R x (R y ) is the extension of the 
QD in the x(y) direction. They obtained the increase of 
J with increasing r for r < 1.3221 and the sharp variation 
of J as a function of interdot detuning for r > S.^The 
influence of an external elec tric field on the exchange 
energy was studied in papers^H for vertically coupled 
self-assembled QD's. Burkard et alP assumed the har- 
monic confinement potential and calculated the exchange 
energy for the vertically coupled QD's using the Hcitlcr- 
London and Hund-MuUiken techniques. These resultsP 
show the monotonic decrease of the exchange energy 
with the increasing in-plane electric field. Pedersen et 
aid calculated the exchange energy for the harmonic 
double-dot confinement potential using the numerically 
exact approach and showed the failure of standard 
approximations (i.e., Heitler-London, Hund-Mulliken, 
Hubbard) even for simple model systems. In the present 
paper, we have elaborated a numerical procedure that 
provides accurate results for the lowest-energy states of 
the two electrons in laterally coupled QD's. Using this 
method, we have performed a systematic study of QD 
nanodevices with different geometry and confinement 
potential profile. In this study, we have applied the 
confinement potentials with a different softness, i.e., we 
have taken into account a variable smoothness of the QD 
interfaced We have investigated a large class of realistic 
potentials with the finite depth: from the soft Gaus sian 
potential to the hard rectangular-like potential! 27 * 29 ! 

A high fidelity of quantum logic operations on spin 
qubits can be achieved if the exchange interaction is pos- 
sibly strong. This leads to the problem of designing such 
a nanodevice, in which the exchange interaction is max- 
imal. A possibility of tuning the exchange interaction 
with the help of external fields is another important is- 
sue of quantum computing in solid-state nanodevices. In 
the present paper, we focus on the electric-field induced 
tuning of the exchange interaction in laterally coupled 
QD's. We note that the elect ronic properties of the elec- 
trostatically gated QD' d 30 | 31 J can be tuned by changing 
the voltages applied to the gates, which changes the po- 
tential confining the electrons in the nanodevice. This 
leads to another method of tuning the exchange interac- 
tion by changing the gate voltages. In the present pa- 
per, we have also investigated this method by studying 
the effect of the variable range and softness of the con- 
finement pot enti al on the exchange energy. In the gate- 
defined QD's,™ the range and softness of the confine- 
ment pot ential ar e determined by the voltages applied to 
the gatesPEUM] 

The present paper is organized as follows: in Section 
2, we briefly describe the theoretical model and the com- 
putational method used in the calculations. Section 3 
contains the numerical results, Section 4 - discussion, 
and Section 5 - conclusions and summary. The details of 
the computational approach are presented in Appendix. 
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FIG. 1: (Color online) (a) Schematic of the coupled QD's. 
Solid (green) curves correspond to the extensions of the QD's. 
(b) Profile of confinement potential energy U plotted as a 
function of x for y = with [solid (red) curve] and without 
[dashed (blue) curve] the external electric field. The vertical 
lines show the boundaries of the left (I) and right (r) electrode. 
The parameters of the nanodevice in panel (b): Ri x = Ri y = 
20 nm, R rx — 20 nm, R ry = 40 nm, d = 80 nm, p = 10, and 
V = 70 mV. 



II. THEORY 

We study the system of two electrons confined in two 
laterally coupled QD's and subjected to a static homo- 
geneous electric field. The lateral QD's are usually cre- 
ated in a quasi-two-dimensional electron gas by apply- 
ing the suitably chosen voltages to the gates, which are 
placed on the surface of the nano device above the plane, 
in which the electrons move HEE1 Therefore, we have as- 
sumed the two-dimensional (2D) motion of the electrons. 
Figure 1 displays (a) the geometry of the nanostructure 
and (b) the confinement potential profile in the x direc- 
tion. We assume that the potential energy of the electron 
in the single QD is described by the power-exponential 
function^ 

U^r) = -U ^e-^- x ^ 2 /<* + (y-y^ 2 / R2 ^ P/2 , (2) 

where index fj, labels the QD's (fi = I and r for the left 
and right QD, respectively), Uo^ is the potential well 
depth (U 0fl > 0), r = (x,y), r Qfl = (x 0fi ,y ^) is the po- 
sition of the QD center, and R^ x (R^y) is the range of 
the confinement potential in the x (y) direction, i.e., it 
determines the extension of the QD in the correspond- 
ing direction. Parameter p (p > 2) describes the softness 
of the confinement potential at the QD boundaries, i.e., 
the smoothness of the QD interfaces.^ For p — 2 we 
deal with the soft Gaussian potential, while for p > 4 
the potential can be treated as "hard", in particular, for 
p — ► oo potential energy ^ takes on the rectangular 
shape. Form ^ of the confinement potential energy al- 
lows us to model a large variety of QD's with a different 
shape, size, and interface smoothness. 

For the coupled QD's the confinement potential is the 
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sum of single QD confinement potentials (|2| 

U conf (r) = U l (r) + U r (r) . (3) 

The x axis is directed along the straight line connecting 
the centers of both the QD's [cf. Fig. 1(a)]. Here, we 
take on xqi — —xq t and yoi — yor = 0. The QD's are sep- 
arated by the potential barrier, i.e., distance d between 
the QD centers is larger than Ri x + R rx . 

We investigate the nanodevice, which consists of the 
coupled QD's placed between the left and right metal 
contacts [cf. Fig. 1(b)]. The electrodes are separated 
by a finite distance L and a static external voltage V 
is applied between them. In the semiconductor region, 
the electrodes generate the homogeneous electric field 
F = (— F, 0, 0), where F = V/L. In the present cal- 
culations, distance L between the boundaries of elec- 
trodes is related to other geometric parameters as follows: 
L > d + 2(Ri x + Rrx)- In electric field F, each electron 
possesses additional potential energy AU(r) given by 

( for x < -L/2 , 

AU(r) = \ -eFx - eV/2 for |z| < L/2 , (4) 
[ -eV for x > L/2 . 

In the present calculations, we measure the energy with 
respect to the electrochemical potential /i/ of the l eft co n- 
tact, i.e., we set fxi = 0. On the contrary to papers i n 
which the infinite range of electric field F is assumed, we 
assume a more realistic space distribution of the electric 
field with the finite range. Formula Q gives the profile of 
the electron potential energy in the electric field created 
by external voltage V, which - in the real nanodevices - is 
applied between the source and drain contacts separated 
by the finite distance. Formulas p]), and Q set up 
a model of the nanodevice (Fig. 1)7 which consists of the 
left (I) and right (r) metal electrodes, and the semicon- 
ductor material, in which both the QD's are embedded. 
The QD's are separated by the barrier potential region. 
The total potential energy of the single electron is given 

by 

U(r)=U conf (r)+AU(r) . (5) 

In the effective mass approximation, the Hamiltonian 
of the two-electron system in the coupled QD's reads 

H=h 1 + h 2 + — — , (6) 

4:TTe e s ri2 

where hj (j = 1,2) is the one-electron Hamiltonian, 
£o is the electric permittivity of the vacuum, e s is the 
static relative electric permittivity of the semiconductor, 
T12 = |r i — r2 1 is the electron-electron distance, and is 
the position vector of the jth electron. The one-electron 
Hamiltonian has the form 

h^-^ + UM, (7) 



where m e is the electron effective band mass. We assume 
that the electron effective mass and the static electric per- 
mittivity do not change across the QD boundaries. This 
assumption is we ll satisfied for the GaAs-based electro- 
static QD's™ 

We solve the two-electron eigenvalue problem by a con- 
figuration interaction (CI) method, which is performed in 
few steps. First, we find one-electron orbital wave func- 
tions <j) v {x, y) using the expansion in a multicenter Gaus- 
sian basis (see Appendix). In the second step, we trans- 
form the one-electron orbitals 4>^(x,y) into the discrete 
representation 0™ n = § v (x m , y n ) on the two-dimensional 
grid (x m ,y n ). More details of this method are given in 
Appendix. Augmenting the one-electron orbitals by the 
eigcnfunctions Xa of the z component of the electron spin 
we obtain one-electron spin-orbitals tp™™ , where v is the 
set of orbital quantum numbers and a is the spin quan- 
tum number. Spin-orbitals ij)™™ are used to a construc- 
tion of Slater determinants. In the final step, we con- 
struct the two-electron wave function as a linear com- 
bination of Ns Slater determinants and solve the two- 
electron eigenvalue equation by the exact diagonalisa- 
tion. All the potential energy matrix elements (including 
the electron-electron interaction energy) have been calcu- 
lated with a high precision by the numerical quadrature 
subroutines. We have performed test calculations for N$ 
= 64, 81, 100, 144, and 169 and obtained a good conver- 
gence for the lowest-energy levels. A good compromise 
between the numerical accuracy and computer time has 
been found for Ns — 81; therefore, the majority of calcu- 
lations has been performed with 81 Slater determinants. 
Finally, we calculate the lowest singlet (Es) and triplet 
(Et) energy levels, and exchange interaction energy J 
[Eq. 0]. In the calculations, we have used the material 
parameters of GaAs, i.e., e s = 12 A and m e = 0.067m e o, 
where m e o is the free electron rest mass, and fix the depth 
(Uqi = Uor = 30 meV) of the confinement potential. The 
present calculations have been performed for circular and 
elliptic QD's with aspect ratios^ r = 0.5, 1, and 2. The 
exchange energy has been calculated as a function of ex- 
ternal electric field F for different shapes, sizes, and ge- 
ometric configurations of the coupled QD's. 



III. RESULTS 

Two lateral QD's form 16 geometric configurations 
that differ from each other by their relative size 
(large/small QD), shape (circular/elliptic QD), position 
with respect to the electrodes (left /right QD), and orien- 
tation with respect to the electric field (ellipse axis par- 
allel/perpendicular to F). In this paper, we present the 
results for the four most characteristic configurations (cf. 
insets of Figs. 2-5). The preliminary results for the two 
circular QD's with the same radius have been presented 
in paperPS Figures 2-4 show the results for the nanode- 
vices, in which the right QD is larger than the left one. 
We remind that the left (right) QD is located near the 
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FIG. 2: (Color online) Exchange energy J as a function of 
electric field F and softness parameter p for two coupled cir- 
cular QD's with Ri x = Ri y = 20 nm, R rx = R ry — 40 nm, 
and d = 80 nm. Inset: schematic of the nanodevice. 



electrode with the higher (lower) potential energy of the 
electron (cf. Fig. 1). These results (Figs. 2-4) have 
been obtained for the left circular QD with fixed size, 
i.e., Ri x = Ri y = 20 nm, and for the different shapes 
and sizes of the right QD: circular (Fig. 2), y-elongated 
elliptic (Fig. 3), and x-elongated elliptic (Fig. 4). We 
have found that - in these nanodevices (cf. the insets of 
Figs. 2-4) - the general electric-field dependence of the 
exchange energy is similar. In the low-field regime, the 
exchange energy takes on either small (Fig. 2) or zero 
(Figs. 3 and 4) values, at higher fields, increases with 
the electric field, and - in the intermediate-field regime - 
exhibits a cusp followed by a broad plateau region (Figs. 
2 and 3), in which J(F) takes on maximal values. At the 
sufficiently high electric field, the J(F) curve possesses 
the second cusp and rapidly falls down to zero. Figures 
2-4 also show another general property of the exchange 
energy: the maximal values J ma x, reached in the plateau 
region, increase with with decreasing p, i.e., with the in- 
creasing softness of the confinement potential. However, 
the detailed J(F) dependence is different for each of the 
nanodevice geometry considered. 

The details of the low-field exchange energy behaviour 
are different in the nanodevices with circular (Fig. 2) and 
elliptic (Figs. 3-5) QD's. For the circular QD's (Fig. 2) 
the exchange energy is non-zero in the absence of electric 
field. In the nanodevice shown in the inset of Fig. 2, the 
electrons in the singlet state occupy the right QD with 
a quite large probability already for F — 0. Even the 
weak electric field causes that both the electrons become 
entirely localized in the right QD, i.e., the double QD 
system starts to act as the single QD. If one of the QD's 
is elliptic (Figs. 3-5), the exchange interaction vanishes 
in the low-field regime, i.e., for < F < F c q, becomes 
non-zero at F — F c q, and increases linearly with F for 
F c0 < F < F cl . At F = F cl the J(F) dependence ex- 
hibits the first cusp, above which J(F) is nearly constant 
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FIG. 3: (Color online) Exchange energy J as a function of 
electric field F and softness parameter p for coupled circu- 
lar and elliptic QD's with Ri x — Ri y = 20 nm, R rx — 20 
nm, R ry = 40 nm, and d = 80 nm. Inset: schematic of the 
nanodevice. 
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FIG. 4: (Color online) Exchange energy J as a function of 
electric field F and softness parameter p for coupled circu- 
lar and elliptic QD's with Ri x — Ri y = 20 nm, R rx — 40 
nm, R ry — 20 nm, and d = 80 nm. Inset: schematic of the 
nanodevice. 

(cf. the plateau regions in Figs. 2 and 3) or changes 
slowly with the electric field (cf. Fig. 4). The exchange 
energy reaches maximal values J m ax for F c \ < F < F C 2, 
exhibits the second cusp at F = F c %, and rapidly van- 
ishes for F > F C 2- For the nanodevice with the larger 
x-elongated QD (Fig. 4) the exchange energy shows the 
variable behaviour in the interval F c \ < F < F C 2 depend- 
ing on the confinement potential softness. According to 
Fig. 4, J decreases with increasing F for p = 2, is nearly 
constant for p — 3, and increases with increasing F for 
p > 4. In the nanodevices shown in the insets of Figs. 3 
and 4, in the low-field regime, the electrons are localized 
in different QD's, i.e., the overlap of the corresponding 
one-electron wave functions vanishes, which leads to the 
vanishing exchange interaction. If the electric field ex- 
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FIG. 5: (Color online) Exchange energy J as a function 
of electric field F and softness parameter p for coupled cir- 
cular and elliptic QD's with Ri x — 20 nm, Ri y — 40 nm, 
nm, and d = 80 nm. Inset: schematic of the 

nanodevice. 
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ceeds the critical value F c q, the electrons in the singlet 
state are entirely localized in the right QD, while the 
electrons in the triplet state become more and more lo- 
calized in the right QD. This leads to the increase of the 
exchange energy in the interval F c q < F < F c ±. The 
plateaus on the J(F) dependence (Figs. 2-4) result from 
the fact that both the electrons are localized in the right 
QD and this localization is almost unchanged in the in- 
terval F c i < F < F C 2- For the sufficiently strong electric 
field the exchange interaction is equal to zero since one 
of the electrons tunnels out of the QD system and is ab- 
sorbed in the right electron reservoir. 

Figure 5 displays the results for the nanodevice with 
the left QD larger than the right one. For F < F c o there 
is no exchange interaction. If the electric field F exceeds 
F c q, the exchange energy becomes non-zero and increases 
as a linear function of the electric field. After reaching 
the maximum at F = F c i, the exchange energy exhibits 
the sharp cusp and falls down to zero. In the nanodevice 
shown in the inset of Fig. 5, the plateau region does not 
exist, which means that F c \ = F c2 . 

The results of Figs. 2-5 show that in the low- and 
intermediate-field regime - the exchange energy is a lin- 
ear function of the electric field and can be parametrized 
as follows: 



J(F) =aF + f3 . 



(8) 



In the nanodevices depicted in the insets of Figs. 3-5, the 
linear parametrization ^ is valid in the electric-field in- 
terval AFn near = F c \ —F c q. Parameter (3 depends on the 
softness of the confinement potential and the geometry of 
the nanodevice. In general, (3 increases with increasing p 
and - for the nanodevice shown in the inset of Fig. 5 - 
takes on the values from —8.56 meV for p = 2 to —6.86 
meV for p = 100. The results of Figs. 3-5 show that 
parameter a is independent of the confinement potential 



FIG. 6: (Color online) Electron density for the triplet state 
[solid (red) curve] and confinement-potential profile [dashed 
(blue) curve] as functions of x coordinate for y = 0. The 
horizontal (green) lines correspond to the three lowest one- 
electron energy levels. Dotted curves show the exchange en- 
ergy J as a function of electric field F (cf. Fig. 3). The 
vertical (orange) arrows show the electric field F for which 
the plots of electron density and confinement potential are 
drawn. In panel (d), the dotted (red) curve depicts the elec- 
tron density of the electron just after the tunneling through 
the triangular barrier. The left QD has the circular shape 
with Ri x = Ri y = 20 nm, and the right QD has the elliptic 
shape with R rx — 20 nm and R ry = 40 nm (cf. inset of Fig. 
3). All the plots are drawn for p = 10 and d — 80 nm. 



softness and the geometry of the nanodevice. It takes on 
nearly constant value a ~ 7.73 [meV/(kV/cm)] for all 
the nanodevices studied in the present work. The physi- 
cal interpretation and possible applications of the linear 
dependence [Eq. (|8|] will be discussed in Section 4. 

The results presented in Figs. 2-5 can be explained 
if we consider the spatial localization of electrons. It is 
convenient to illustrate the distribution of the electrons 
in the coupled QD's with the help of the electron density 
f>(r) defined as 



e(r)=£)<*|5(r-r j )|*) 



(9) 



where 'F = ^(ri,^) is the two-electron wave function. 

Solid (red) curves in Figs. 6 and 7 display the cross sec- 
tions of the electron density g(x, y) by the plane y = for 
the triplet state. The corresponding profiles of the one- 
electron potential energy are plotted by dashed (blue) 
curves in Figs. 6 and 7. The horizontal (green) lines show 
the three lowest-energy levels E v of one-electron states, 
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FIG. 7: (Color online) One-electron density for the triplet 
state [solid (red) curve] and confinement-potential profile 
[dashed (blue) curve] as functions of x coordinate for y — 0. 
The horizontal (green) lines correspond to the three lowest 
one-electron energy levels. Dotted curves display the ex- 
change energy J as a function of electric field F (cf. Fig. 5). 
The vertical (orange) arrows show the electric field for which 
the plots of probability density and confinement potential are 
drawn. In panel (d), the dotted (red) curve depicts the prob- 
ability density of the electron just after the tunneling through 
the triangular barrier. The left QD has the elliptic shape with 
Ri x = 20 nm and Ri y = 40 nm, the right QD has the circular 
shape with R rx — R ry — 20 nm (cf. inset of Fig. 5). All the 
plots are drawn for p = 10 and d = 80 nm. 



which mostly contribute to the two-electron triplet wave 
function. In Figs. 6 and 7, we have also presented the 
J(F) dependencies obtained for these nanodevices (cf. 
Figs. 3 and 5 for p — 10). The vertical (orange) arrows 
show the values of the electric field, for which the plots 
of g, U, and E v are drawn. Figures 8 and 9 display the 
contours of the electron density on the x — y plane for the 
singlet and triplet states for the same two nanodevices. 
The results presented in Figs. 6 and 8 correspond to the 
nanodevice, in which the y extension of the right QD is 
two times larger than that of the left QD. We see that 
the smaller values of g(x, 0) in the larger (right) QD [cf. 
Figs. 6(a,b)] are compensated by the larger spreading of 
the electrons in the y direction (cf. the right panel of Fig. 
8). 

The results displayed in Figs. 6-9 allow us to trace the 
changes of electron localization in the two-electron sys- 
tem confined in the two coupled QD's, which result from 
the action of the external electric field. In the absence 
of the electric field, the electrons are localized in the dif- 
ferent QD's and there is no overlap between their wave 
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FIG. 8: Contours of electron density for singlet (left panels) 
and triplet (right panels) states on the x — y plane for different 
values of the external electric field F (in kV/cm). The darker 
color corresponds to the larger electron density. Thin solid 
curves correspond to the sizes of the QD's. The left QD has 
the circular shape with Ri x — Ri y = 20 nm and the right QD 



has the elliptic shape with R r 



20 nm, R r 



40 nm; the 



other parameters take on the values p = 10 and d = 80 nm. 



functions [cf. Figs. 6(a) and 7(a), and Figs. 8 and 9 
for F = 0]. Therefore, the exchange interaction between 
the electrons vanishes. If we apply the external electric 
field, the electrons start to tunnel through the potential 
barrier from the left to the right QD and interact via the 
exchange coupling [Figs. 6(b) and 7(b)]. In this regime, 
the increasing electric field causes the fast linear increase 
of the exchange energy [Figs. 6(b) and 7(b)]. We note 
that - even in the low electric field regime the electron 
distribution considerably changes in the singlet state, but 
is only slightly distorted in the triplet state (cf. Figs. 8 
and 9 for F = and F > 0). If the right QD is suffi- 
ciently large, the increasing electric field causes that both 
the electrons become localized in this QD in cither spin 
state [cf. Fig. 6(c) and Fig. 8 for F — 1.52 kV/cm]. In 
the singlet state, the electrons are localized in the central 
region of the right QD, while in the triplet state, the elec- 
tron density exhibits two maxima clearly separated in the 
y direction (cf. Fig. 8 for F = 1.52 kV/cm). In this field 
regime, the exchange energy reaches the largest values. 
The double occupancy of the right QD and the separation 
of the triplet electrons in the y-direction do not change 
in a rather broad range of the electric field, which leads 
to the broad plateau region (Fig. 3). If the electric field 
exceeds the critical value F C 2, one of the electrons tun- 
nels through the right triangular barrier out of the right 
QD and the exchange interaction rapidly falls down to 
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FIG. 9: Contours of electron density for singlet (left panels) 
and triplet (right panels) states on the x — y plane for different 
values of the external electric field F (in kV/cm). The darker 
color corresponds to the larger electron density. Thin solid 
curves correspond to the sizes of the QD's. The left QD has 
the elliptic shape with Ri x = 20 nm and Ri y = 40 nm, the 
right QD has the circular shape with R rx — R ry — 20 nm; 
the other parameters take on the values p — 10 and d — 80 
nm. 



zero [cf. Figs. 6 (d) and 8 for F — 3.16 kV/cm and 
Figs. 7(d) and 9 for F = 2.40 kV/cm]. The dotted (red) 
curves in the right part of Figs. 6(d) and 7(d) illustrate 
the electron density just after the tunneling. The tunnel- 
ing electron is absorbed in the electron reservoir of the 
right electrode. Let us note the corresponding changes of 
the one-electron energy levels (cf. the horizontal lines in 
Figs. 6 and 7). For F < F c2 the one-electron energy lev- 
els lie below the right-electrode continuum energy edge, 
i.e., electrochemical potential \i r = — eV of the right elec- 
trode. This means that these electron states are bound. 
For F — F c2 the one-electron energy levels reach \i r , i.e., 
the electrons cease to be bound and form resonant states. 
Therefore, in the strong-field regime, the electrons tun- 
nel via these resonant states through the right triangular 
barrier [Figs. 6(d) and 7(d)] to the right contact. 

The nanodevice with the left QD larger than the right 
one (Figs. 7 and 9) shows the similar electric-field be- 
haviour to that obtained for the nanodevice with the left 
QD smaller than the right one (Figs. 6 and 8) in the 
regime of low and intermediate electric fields. We re- 
mind that - in this regime - the exchange energy is zero 
for F < F cQ and grows linearly with F for F c0 < F < F cl . 
However, after reaching the maximum value J max for 
F = F c i, the exchange energy rapidly falls down to zero, 
i.e., the J(F) dependence becomes qualitatively different 
from that shown in Figs. 3 and 6. This behaviour results 
from the rapid change in electron localization that occurs 



at F = F c i ■ For F < F c i the electrons in the triplet state 
are localized in the different QD's [cf. Figs. 7(b,c)]. For 
the triplet state, the increase of the electric field above 
F c i does not generate the double occupancy of the right 
QD, like in the nanodevice shown in Fig. 3, but leads to 
the immediate tunneling of one of electrons to the right 
electron reservoir [Fig. 7(d)]. In the electric-field regime 
F > F c i, the right QD in this nanodevice can not be 
occupied by the two electrons in the triplet state, which 
causes that the exchange interaction vanishes. Fig. 7(d) 
shows that for the sufficiently strong electric field the en- 
ergy of the first excited one-electron state exceeds the 
electrochemical potential of the right contact, i.e., the 
resonant state is formed. In this case, we are dealing with 
the resonant tunneling via the first excited one-electron 
state. We note that this one-electron state yields the 
large contribution to the triplet two-electron wave func- 
tion. 

We have also studied the dependence of the exchange 
energy on the softness of the confinement potential. Fig- 
ure 10 shows maximum exchange energy J max as a func- 
tion of softness parameter p. The maximal value J max is 
taken for F = 0.7576 kV/cm, i.e., in the plateau region, 
for the nanodevice shown in the inset of Fig. 3. We see 
that Jmax decreases if p increases, i.e., if the confinement 
potential becomes more hard. This dependence can be 
approximated by the exponential function 

J m ax(p) = Ai exp(-Cip) + Si , (10) 

where Ai = 2.4740 meV, B x = 1.4059 meV, and C x 
I i. 2(H) 1 . The exponential parametrization [Eq. ( |To| )] re- 
sults from the fact that - for the fixed confinement- 
potential ranges - the effective quantum capacity of the 
QD's increases with increasing p (cf. inset of Fig. 10). 
If p increases, the electrons localized in the right QD be- 
come more separated from each other, which leads to the 
exponential decrease of the overlap of electron wave func- 
tions, which in turn gives rise to the exponential decrease 
of the exchange energy. 

For a possible experimental realization of the model 
nanodevices studied in the present paper it is interesting 
to find a direct dependence of the exchange energy on the 
size of the nanodevice. For this purpose we have calcu- 
lated the maximum exchange energy when scaling all the 
linear dimensions of the coupled QD system. We consider 
the nanodevice displayed in the inset of Fig. 3, for which 
the exchange energy takes on the maximal values in the 
broad plateau region (cf. Fig. 3 for p= 10). We have de- 
fined the size scaling factor as s = R ac tuai / RiniUah i-e., s 
is equal to the ratio of the actual linear dimension Ractuai 
to its initial value Rmitiai- As the reference nanodevice 
with the initial values of the linear dimensions we take 
that with Ri x = Ri y — 20 nm, R rx = 20 nm, R ry = 40 
nm, and d = 80 nm (cf. Fig. 3). The calculations of J m ax 
have been performed for the set of nanodevices charac- 
terized by s times enlarged confinement potential ranges, 
i.e., sRi xl sRi yi sR rx , and sR ry and interdot distance sd. 
In the calculations, we fix the strength F of the electric 
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1.7874 (cf. solid curve in Fig. 11). The exponential de- 
pendence [Eq. (Ill] is similar to that given by Eq. (10 1 



FIG. 10: (Color online) Maximum exchange energy Jmax as 
a function of softness parameter p for electric field F = 0.7576 
kV/cm. Dots show the results of numerical calculations 
and the solid curve shows the fitted exponential function 
[Eq. (| 10p ] . The parameters of the nanodevice: Ri x = Ri y = 20 
nm, Rrx = 20 nm, R ry = 40 nm, and d = 80 nm. Inset: Total 
potential energy U of the electron as a function of x and p for 
y = and for fixed F = 0.7576 kV/cm.. 





-500 500 1000 

x [nm] 



F = 0.7576 kV/cm 



FIG. 11: (Color online) Maximum exchange energy J ma x as 
a function of size scaling factor s for p = 10 and for electric 
field F — 0.7576 kV/cm. Dots show the results of numerical 
calculations and the solid curve shows the fitted exponential 
function [Eq. (111]. The initial sizes (s = 1) of the QD nan- 
odevice are Ri x = Ri y = 20 nm, R rx = 20 nm, R ry — 40 
nm, and d = 80 nm. Inset: Total potential energy U of the 
electron as a function of x and s for y = and for fixed F = 
0.7576 kV/cm. 



field, i.e., we have to scale accordingly the interelectrode 
distance L — > sL and the applied voltage V — ■*■ sV. The 
numerical results are displayed in Fig. 11 by the full 
(red) dots. In the interval 1 < s < 5 these results can be 
parametrized by the exponential function 

J max {s) = J 4 2 exp(-C 2 s) + B 2 , (11) 

where A 2 = 9.9040 meV, B 2 = 0.0710 meV, and C 2 = 



and can also be interpreted as resulting from the size 
effect. If the total size of the nanodevice increases (cf. 
inset of Fig. 11), the overlap between the one-electron 
wave function decreases exponentially with increasing s. 
Also the localization of electrons in the QD's becomes 
weaker if the total size of the nanodevice grows. We note 
that parametrization (11) is valid for s < 5 only. If the 
size of the nanodevice is sufficiently large, i.e., the size 
scaling factor exceeds ~ 5, the exchange energy rapidly 
falls down to zero. The disappearance of the exchange 
interaction in the large-size nanodevice results from the 
derealization of electrons, which can be explained using 
the potential energy profiles (cf. inset of Fig. 11). In 
order to keep the electric field constant when enlarging 
the nanodevice size s times we have to apply the s times 
higher voltage. This leads to the lowering of electrochem- 
ical potential p r = —eV of the right contact. Simulta- 
neously, the energy of the electrons localized in the right 
QD grows with respect to fi r . We have checked that for 
s ~ 5 the first excited-state one-electron energy level ex- 
ceeds the electrochemical potential of the right contact. 
Moreover, the triangular barrier near the right contact 
becomes more and more penetrable for the electrons if 
the size of the nanodevice increases. In these conditions, 
one of the electron tunnels out of the right QD to the 
right reservoir and the exchange interaction vanishes. 



IV. DISCUSSION 

The results of Figs. 3-5 show that - in the nanode- 
vice, which consists of the elliptic QD - the static homo- 
geneous electric field applied in the coupled QD's region 
can switch on and off the exchange interaction. In the 
nanodevice, which consists of the two circular QD's (Fig. 

2) , the exchange interaction is non-zero at F = and the 
electric field can only switch off the exchange interaction. 
This behaviour (Fig. 2) is similar to that observed in the 
single QDp^ In the nanodevices shown in the insets of 
Figs. 2-4, the exchange energy exhibits a plateau in a 
broad electric-field regime. This plateau ends up at crit- 
ical electric field F c2 , above which the exchange energy 
rapidly vanishes. If the electric field exceeds F c2 , one of 
the triplet electrons tunnels from the right QD to the 
right contact and the triplet state ceases to be bound. 
Therefore, in this field regime, we can not speak about 
the exchange interaction. The electrons in the singlet 
state become unbound if the electric field exceeds F c2 by 
an amount AF$- We have found that AF$ — AFu near , 
i.e., AFs is approximately equal to the width of the linear 
J(F) dependence [Eq. 

The critical electric field F c2 increases with increas- 
ing p (cf. Figs. 2-4). Simultaneously, the increasing p 
leads to the decreasing maximum value of the exchange 
energy reached in the plateau region (cf. Figs. 2 and 

3) . Both these effects result from the increasing effective 



9 



size of the QD's. For fixed parameters Ri x ,Ri y , R rx , and 
R ry the effective size of the QD increases with increasing 
p, i.e., increasing hardness of the confinement potential 
(cf. inset of Fig. 10). This leads to the decreasing over- 
lap between the electron wave functions and the weaker 
electron localization, which in turn causes the decline of 
the exchange energy. Moreover, if the effective size of 
the QD's is larger, we have to apply the stronger electric 
field in order to liberate one of the electrons from the 
right QD , which gives rise to the increase of F C 2 . 

The critical electric field F c0 , below which J = and 
above which J > 0, decreases with increasing p (cf. Figs. 
2-5). This dependence results from the decreasing ef- 
fective thickness of the potential barrier separating both 
the QD's with increasing p (cf. inset of Fig. 10). If the 
barrier is thinner, the electrons tunnel through it with 
the larger probability and the right QD becomes doubly 
occupied at the lower electric field. 

In the low- and intermediate-field regime, the nanode- 
vices with the laterally coupled QD's possess an impor- 
tant property: the exchange energy is a linear function of 
the electric field, i.e., it can be conveniently tuned by ap- 
plying the external voltage. In the nanodevice with the 
small right QD (cf. Fig. 5), the increasing electric field 
switches on the exchange interaction at F = F c q, leads 
to the linear increase of J in a broad regime of F, and 
switches it off at F — F c \. Recently, the linear depen- 
dence of the exchange interaction energy on the electric 
field has been found in vertically coupled self-assembled 
QD's. 2b In the present paper, we have obtained this lin- 
ear dependence for the laterally coupled QD's with dif- 
ferent geometric configurations and different softness of 
the confinement potential (cf. Figs. 2-5). Electric-field 
regime AFi inear of this linear dependence is considerably 
broader in the nanodevices that contain the small right 
QD (Fig. 5) than in the nanodevices with the large right 
QD (Figs. 2-4). According to Fig. 5, AFu near extends to 
~ 0.8 kV/cm. For comparison, in Figs. 2-4, AFi inear ~ 
0.2 kV/cm. The larger width of interval AFu near ob- 
tained for the nanodevice depicted in the inset of Fig. 5 
causes that the maximum value J max of the exchange en- 
ergy is considerably larger than those for the nanodevices 
shown in Figs. 2-4. In the former case, J m ax reaches 6 
meV. 

In order to get a more deep physical insight into the lin- 
ear J(F) dependence, we have investigated the behaviour 
of the lowest singlet and triplet energy levels. We have 
found that - in the regime < F < F C 2 - the field de- 
pendence of these energy levels can be very accurately 
parametrized by the linear functions (Fig. 12) 

E S , T (F) = a s , T F + (3 S ,T ■ (12) 

Parameters as, olt, Ps, an d Pt take on different val- 
ues in the different regimes of the electric field, i.e., 
[0,F c oUF c o,F c i],[F c i,F c2 ], that we label (0), (1), and 
(2), respectively, but are independent of F within each 
regime. In the low-field regime, i.e., for < F < F c q, 
— a r an d Ps ~ Pt^ which leads to the ze- 
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FIG. 12: (Color online) Energy of singlet (Esw, Es s ) and 
triplet (Etw,Ets) states and exchange energy J (dash-dotted 
curves) as functions of electric field F for two distances d 
between the QD centers. Inset shows the geometry of the 
nanodevice. Subscripts w and s correspond to the weak (d = 
80 nm) and strong (d = 40 nm) tunnel coupling, respectively. 
The other parameters of the QD's are the same as in Fig. 3 
with p — 10. For clarity the curves with label s are shifted 
upwards by 5 meV. Thin vertical lines show critical fields 
F c q and Fd for the weak tunnel coupling, the vertical arrow 
shows the cusp on the Ets(F) plot, which gives rise to the 
corresponding cusp on the J(F) plot for the strong tunnel 
coupling. Labels (0), (1), and (2) denote the electric-field 
regimes [0,-F c o], [F c o,F cl ], and [F cl ,F c2 ], respectively, for the 
weak tunnel coupling. 



roing of the exchange energy, i.e., the singlet-triplet 
degeneracy. In the intermediate-field regime, i.e., for 
F c o < F < F c i, the singlet-triplet degeneracy is lifted. In 
this field regime, a^ 1 ' < a^ and /3g > . Therefore, 

we obtain a = ct^ — a^ > and j3 = P^ — Pg_ < 0, 
which leads to the linear J(F) dependence [Eq. pj]. 

It is interesting that - for the nanodevices shown in the 
insets of Figs. 2-3 - linear parametrization (l2"| is also 
valid in the regime of rather strong electric fields, i.e., for 

(2) (2) 

F c \ < F < F c2 . The values of parameters a S T and j3 S T 
arc different from those obtained for the lower fields, but 
are approximately constant within this field regime. In 
field regime (2), a^' ~ atp and (3^p > P^, which gives 
rise to the plateau of the exchange energy (cf. Figs. 2-4) 

with J m ax — Pt ~ Ps = cons t- 

We have also found another interesting feature of the 
linear parametrization dll. Parameter a = AJ/AF that 
determines the rate of changes of the exchange energy 
with the electric field in the linear regime takes on almost 
the same values for all the nanodevices described by the 
parameters quoted in the captions of Figs. 2-5. Param- 
eter a is independent of the softness of the confinement 
potential and the geometry of the nanodevice (cf. Figs. 
2-5). Considering all the J(F) dependencies, displayed in 
Figs. 2-5, i.e., studying several different nanodevices, we 
have found that a = 7.73 ± 0.13 [meV/(kV/cm)]. In the 



10 



lateral QD's, the linear J(F) dependence is a non-trivial 
property that can not be explained by the non-degenerate 
first order perturbation theoryPSl This effect occurs in 
the intermediate-field regime, in which the electron wave 
functions are considerably distorted with respect to those 
for F — (cf. Figs. 6-9). In particular, we note the 
rapid change of the singlet-state localization in the low- 
and intermediate-field regime (Figs. 8 and 9). Moreover, 
we have found that even for the same geometric configu- 
ration of the QD's and in the same electric-field regime, 
the linear J(F) dependence disappears if the tunnel cou- 
pling between the QD's is sufficiently strong. In order to 
show this effect, we have considered the two nanodevices 
with the geometric configuration shown in the insets of 
Figs. 3 and 12. Fig. 12 displays the results for the two 
nanodevices characterized by different separations d be- 
tween the QD centers and the same values of all other 
parameters. For d — 80 nm, i.e., for the weak interdot 
tunnel coupling, the singlet and triplet energies as well 
as the exchange energy are linear functions of the elec- 
tric field. We also observe the cusp on the curve Et w (F), 
which leads to the corresponding cusp on the curve J(F) . 
In the low electric-field regime and for the weak tunnel 
coupling, the exchange interaction vanishes due to the 
singlet-triplet degeneracy. For d — 40 nm, i.e., for the 
strong tunnel coupling, a completely different J(F) be- 
haviour has been obtained. In this case, the singlet-state 
energy is a non-linear function of the electric field, which 
leads to the non-linear J(F) dependence. However, the 
triplet energy Ef s is a piece-wise linear function of F 
with the cusp shown by the arrow in Fig. 12. In the low 
electric-field regime and for the strong tunnel coupling, 
the singlet-triplet degeneracy is lifted and the exchange 
interaction is non-zero even for F = 0. 

The linear/non- linear J(F) dependence can be further 
discussed if we consider the changes of the charge grav- 
ity center position, i.e., x c = (xi + x^jl. We have 
calculated the expectation value (x c ) for the weak and 
strong interdot tunnel coupling using the field-dependent 
two-electron wave functions ^(ri,^) for the singlet and 
triplet states. According to Eq. Q, the value 2(x c ) de- 
termines the electric-field contribution to the potential 
energy of two electrons. Figure 13 shows that (x c ) is 
a linear function of F for the triplet states in both the 
cases of the weak and strong tunnel coupling and for the 
singlet state only in the case of the weak tunnel coupling. 
However, for the singlet state and the strong tunnel cou- 
pling the dependence of (x c )s s on the electric field is 
nonlinear. We also observe that the expectation values 
{x c )sw, {x c )tw, and (x c )t s exhibit jumps at certain val- 
ues of the electric field. Besides they are linear functions 
of F. Figure 13 shows that the charge gravity center 
position follows the electron distribution shown in Fig. 
8 for d = 80 nm. At low electric fields, the center of 
charge gravity is localized near the center of the nanode- 
vice, i.e., at x ~ 0, for the triplet as well as the singlet 
states, which results from the single occupancy of both 
the QD's. The increasing electric field only slightly shifts 
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FIG. 13: (Color online) Expectation value (x c ) of the charge 
gravity center position as a function of electric field F for 
the singlet (subscript S, curves with full circles and squares) 
and triplet states (subscript T, curves with open circles and 
squares) for d = 40 nm (subscript s, red curves) and d = 80 
nm (subscript w, blue curves). The nanodevice is the same 
as in Fig. 12. 



(x c ) to the right QD. The single QD occupancy remains 
unchanged for the triplet state up to F = F c \. How- 
ever, for the singlet state the second electron tunnels to 
the right QD at F = F c q, which leads to the jump of 
(x c )sw (Fig. 13). For F > F c \ and the weak tunnel cou- 
pling the charge gravity center is approximately localized 
near the center of the right QD in both the spin states 
(cf. the plots of (x c )sw and (x c }t w in Fig. 13). For the 
weak tunnel coupling, the localization of the electrons in 
the singlet state remains almost unchanged in the regime 
F c o < F < F c ±. For the strong tunnel coupling (x c )s s in 
the singlet state increases as a non-linear function of F 
starting from a non-zero value for F = 0. In the triplet 
state, the jump of (x c )t s at F = F c % originates from the 
rapid change of the localization from the single to double 
occupancy of the right QD. 

The results of Fig. 13 can be translated into the energy 
dependencies displayed in Fig. 12. The linear J(F) de- 
pendence occurs if the interdot tunnel coupling is weak. 
Then, in the triplet state, the single-electron occupancy 
of each QD occurs in electric-field regime [0, F c i\. In the 
field regime [F c0 ,F c i], the right QD is doubly occupied 
by the singlet electrons. If the interdot separation is suffi- 
ciently large, the increasing electric field slowly shifts the 
center of the charge gravity towards the right QD. Si- 
multaneously, the electrons are compressed at the right 
potential barriers of QD's, which hampers this shift. The 
net shift is proportional to the electric field, which results 
in the linear dependencies Es(F) and Et(F). If the in- 
terdot tunnel coupling is sufficiently strong, the electrons 
can tunnel through the barrier with a rather large proba- 
bility, which results from the considerable overlap of the 
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wave functions centered in both the QD's. Therefore, for 
the strong tunnel coupling the increasing electric field 
leads to the non-linear increase of the exchange energy. 



V. CONCLUSIONS AND SUMMARY 

The results of the present paper allow us to discuss the 
effect of external electric fields of different origin on the 
electronic properties of QD nanodevices. The electric 
fields can be created by the different electrodes, which 
surround the QD region. In a direct way, we have in- 
vestigated the effect of the static homogeneous electric 
field, which is usually created by the source and drain 
electrodes. In an indirect way, we have also studied the 
effect of the gate electrodes. We are able to determine the 
effect of the g ates, sin ce the gate-controlled QD's (elec- 
trostatic QD's)^ 31 * 32 lare induced by the inhomogeneous 
electric field, which is created by the gates. Therefore, 
the QD confinement potential, in particular, its shape, 
range, potential-well depth, and softness, are determined 
by the voltages applied to the gates. This means that 
studying the dependence of the electron states on these 
parameters of the confinement potential we are investi- 
gating the effect of the gates. 

In the present paper, we focus on the exchange inter- 
action, which pl ays an important role in a manipulation 
with spin qubitsP^lMI The exchange energy can be effec- 
tively tuned by changing the external electric field, i.e., 
changing the voltages applied to the electrodes. By in- 
creasing the homogeneous electric field (bias voltage) we 
can switch on/off the exchange interaction. The critical 
values of the electric field, for which this on/off switch- 
ing occurs, depend on the softness of the confinement 
potential, i.e., can be changed by changing the gate volt- 
ages. We have found the existence of the plateau in the 
exchange energy versus electric field dependence. This 
plateau occurs if the QD, into which the electrons are 
pushed by the electric field, is larger than the other QD. 
In the plateau regime, the exchange energy take on is 
maximal values. Therefore, the on/off switching the ex- 
change interaction occurs between zero and maximum 
value. We have determined the critical electric fields and 
optimal nanodevice parameters, for which this switching 
is the most effective. 

We have shown that - in the nanodevices with the 
weakly coupled lateral QD's at moderate electric fields 
- the exchange energy is a linear function of the electric 
field. We have found that the parameter a = AJ/AF 
that determines the rate of changes of J(F) is nearly 
the same for different nanodevices. The constancy of 
a suggests that this parameter is universal for a large 
class of nanodevices based on the laterally coupled QD 's 
provided that the interdot tunnel coupling is weak. We 
have also demonstrated that for the sufficiently strong 
interdot tunnel coupling the J(F) dependence becomes 
non-linear. 

The numerical approach proposed (see Appendix) is 



convenient for the calculations of few-electron states in 
QD's in the external electric field. Among several advan- 
tages (cf. Appendix) of this approach, we would like to 
underly the most important one. Namely, this numeri- 
cal procedure possesses the following physical property: 
in the low- and intermediate- field regime, it allows us to 
describe bound states of electrons in the electric field of 
finite range, which acts in the real nanodevices. Let us 
mention that in the electric field of in finite range, usually 
assumed in the papers on this subject J2U2SI we always deal 
with non-stationary states that - for weak fields only - 
can be treated as quasi-bound. 

Having at disposal the nanodevices based on the gate- 
controlled lateral QD's, we can intentionally change the 
shape, range, and potential- well depth of the confinement 
potential by changing the voltages applied to the gates. 
The manipulation of electronic states in laterally coupled 
gate-controlled QD's, performed by changing the gate 
voltages, is very effective. Therefore, the nanodevices 
based on the lateral QD's are very promising in quantum 
computing based on spin qubits. In the present paper, 
we have shown how to control and tune the exchange 
interaction by changing the external electric field. Our 
results should be helpful in designing the nanodevices for 
the spin qubit processing with the help of the controlled 
exchange interaction. 
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APPENDIX: COMPUTATIONAL METHOD 

We give a brief description of the numerical method 
proposed in the present paper in order to obtain the 
one-electron wave functions in the discrete representation 
that is convenient in the two-electron calculations. First, 
we set up the computational box for given geometric pa- 
rameters of the QD nanodevice. Usually, we take on its 
x extension L x = 400 nm and y extension L y — 140 nm. 
Next, we solve the one-electron Schrodinger equation by 
variational means. For this purpose we apply the multi- 
center Gaussian basis defined inside the computational 
box 



9i 3 {x,y) = exp[-7x(x - x^ 2 - -y 2 (y ~ Vjf 



(A.l) 



where 71 and 72 are the variational parameters. 

The centers of Gaussians (xi,yj), where i — I,..., I 
and j — 1, . . . , J, form a grid within a large rectangle, 
which encompasses both the QD's. i.e., the edges of the 
rectangle are longer than the double range of the confine- 
ment potential in each direction (cf. Fig. 14) . Gaussians 
(A.l I are chosen to cover the region of electron localiza- 



tion in the QD's and the region of the right electrode, 
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FIG. 14: (Color online) (a) Profile of electron potential en- 
ergy U as a function of x for y = and F = 1.2625 kV/cm. 
Symbols I and r denote the left and right electrode, respec- 
tively. Vertical dashed lines correspond to the boundaries 
of the electrodes, (b) Schematic of the computational box. 
Thick (red) dots display the positions of the centers of Gaus- 
sians and thin (black) dots show the grid points used in dis- 
crete representation (A.3I. The nanodevice parameters used 
in the plots: Ri x = Ri y = 20 nm, R rx = 40 nm, R ry — 40 
nm, p = 10, and d = 80 nm. 



to which the electrons are shifted. The variational wave 
function <f> v {_x,y) for one-electron state v is taken on in 
the form of a linear combination of Gaussians ( A.l I 



/ J 

b v (x, y) = C v '•;',,</<.;':•'•• y) 

i=l j = l 



(A.2) 



where C v is the normalization constant. The values of 
I and J are chosen according to the actual size of the 
nanodevice and extend up to I max — 50 and J max — 33. 
This means that we have at disposal 1650 basis elements 
(A.l I, which assures the high accuracy of one-electron 
solutions. The matrix elements of one-electron Hamil- 



tonian ^ are calculated in basis (A.l I as follows: the 
matrix elements of the kinetic energy are calculated an- 
alytically and the matrix elements of the potential en- 
ergy are calculated by a numerical quadrature. Solving 
the generalized eigenvalue problem of Hamiltonian ^ 
in basis (A.l I we obtain linear parameters and energy 
eigenvalues k v for one-electron states v = 1, . . . , N v . The 
values of non-linear variational parameters 71 and 72 are 
determined from the minimization of the ground-state 
energy Eq. We have checked that 71 and 72 change only 
slightly when minimizing the excited-state energy levels; 
therefore, we take 71 and 72 to be the same for each state 
v. 

In the last step, we define the fine grid (x m , y n ), where 
m = 1, . . . , M and n = 1, . . . , N, and find the discrete 
representation of the one-electron wave functions by set- 
ting the M x N matrix 



(A.3) 



The one-electron wave functions in representation (A.3 1 



are used to construct the Slater determinants and next to 
perform the CI calculations. In the two-electron calcula- 
tions, we have used N u = 13 one-electron states to con- 
struct up to Ng = 169 Slater determinants. The number 
of mesh points was M x N = 123 x 91. 

We would like to emphasize the following advantages of 
the present approach: (i) It allows us to solve the electron 
eigenproblem in the real nanodevice, in which the elec- 
tric field has the finite range, i.e., the voltage is applied 
between the electrodes separated by the finite distance. 
This numerical feature possesses an important physical 
consequence: if the electric field is not too strong, we are 
dealing with the well-defined stationary bound states. It 
is commonly assumed that the electric field possesses the 
infinite range, i.e., instead of the form Q the potential 
energy has the form AU — —eFx for — 00 < x < +00. In 
this case, the electron states are unbound for arbitrary 
F. This assumption leads to serious problems in the 
physical interpretation and does not allow us to describe 
rigorously the real finite-size nanodevices. (ii) It can be 
applied to arbitrary confinement potential and inhomo- 
geneous electric field. In particular, this approach can be 
easily extended to multiple coupled QD's with more than 
two QD's. (iii) The present approach can be extended to 
few-electron systems with more than two electrons and 
to three-dimensional (3D) nanostructures. 
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